Magnetism and local distortions near carbon impurity in 7-iron 
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Local perturbations of crystal and magnetic structure of 7-iron near carbon interstitial impurity 
is investigated by ab initio electronic structure calculations. It is shown that the carbon impurity 
creates locally a region of ferromagnetic (FM) ordering with substantial tetragonal distortions. 
Exchange integrals and solution enthalpy are calculated, the latter being in a very good agreement 
with experimental data. Effect of the local distortions on the carbon-carbon interactions in 7-iron 
is discussed. 

PACS numbers: 61.72.Ji, 75.20.En, 75.30.Et, 71.20.Be 
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Steel is a material playing unique role in our civiliza- 
tion. Its main chemical composition is very simple, this 
is just iron and carbon. Nevertheless, surprisingly, ba- 
sic microscopic physics determining phase and structural 
states of the steel is rather poorly understood yet. In par- 
ticular, mechanisms of development of lattice instability 
and martensitic transformation at cooling down of 7-Fe 
are still unknown. It is commonly accepted now that 
magnetism is of crucial importance for the phase stabil- 
ity of iron and its alloys |3, 0], however, more or less 
detailed theoretical studies have been carried out only 
for a particular case of Invar Fe-Ni alloys [3, |4|, [5j • 

Commonly used steel is based on the low-temperature 
(a) bec phase of iron. However, morphology of its mi- 
crostructure which is decisive for all practical applica- 
tions is formed during the quenching process from high- 
temperature fee (7) phase. Kinetics of the 7 — 0: transi- 
tion is very sensitive to carbon concentration. The state 
of carbon in a-iron was a subject of numerous investi- 
gations @, 0, HI • The state of carbon in 7-iron is much 
less understood, and even its solution enthalpy calculated 
by state-of-art ab initio approach strongly disagrees with 
the experimental data Q which is quite unusual. 

The most probably, magnetic state of 7-iron is rele- 
vant for structural properties of its alloys. Its magnetic 
state is strongly frustrated which leads to existence of nu- 
merous c omp licated magnetic structures with very close 



energies (lfj, |12|. Role of the lability of magnetic 



structure and the frustrations is discussed already in the 
context of Invar problem 0, [l| . Here we present the re- 
sults of ab initio calculations of the electronic structure, 
lattice and magnetic properties of carbon solid solution 



in 7-iron. It turns out that the carbon interstitial in oc- 
tahedron void results in an essential local magnetic po- 
larization and strong lattice distortions which should be 
taken into account, in particular, to obtain correct value 
of the solution enthalpy. 

We used the SIESTA package of first-principles elec- 
tronic structure calculations [13l . Il4| with the generalized 
gradient approximation for the density functional [la ]. 
Earlier the same approach has been successfully used to 
calculate various properties of bulk and surface iron [lH ] 
as well as Fe clusters 17| . To calculate exchange inter- 
actions of the effective Heisenberg model 



H, 



eff 



i,3 



(Si are the classical spins defined by the direction and 
magnitude of obtained magnetic moments) a standard 
density functional approach has been used based on the 
"magnetic force theorem" [3]. We optimized first the 
structure and then use the implementation of the Green's 
function [H,^ into the LMTO method [1H to calculate 
the effective exchange parameters which is not possible in 
the framework of SIESTA. Justification of this approach 
will be presented below, at the discussion of the Table HT1 
To check accuracy of our approach we first studied 
structural relaxation effects in pure 7-Fe which is known 
to result in tetragonal deformations of initial fee lat- 
The computational results presented in 



ticc 



22, 22 



the Table U demonstrate a reasonable agreement with 
the previous calculations and with experimental data for 
thin films of 7-Fe 22| . We have considered the following 
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TABLE I: 

Lattice parameters, magnetic moments, and energies of differ- 
ent magnetic configurations per atom for 7-Fe with atomic re- 
laxation taken into account; numbers in parenthesis are taken 
from ReflU 





FM 


AFM 


AFMD 


a, A 


3.58 (3.45) 


3.44 (3.45) 


3.57 (3.49) 


c/a 


1.08 (1.18) 


1.09 (1.09) 


1.05 (1.09) 


M, fjL B 


2.5 


1.8 


2.3 


E - Eafmd, meV 


33 


45 






magnetic structures: ferromagnetic (FM), antiferromag- 
netic with the staggered magnetization in < 001 > direc- 
tion (AFM), and double antiferromagnetic (AFMD), or 

"+H " (see Fig. [IJ . The latter magnetic configuration 

is one of most energetically favorable for 7-Fe pU. H^. 




FIG. 1: (color online) Fragment of crystal and magnetic struc- 
ture of 7-Fe for the AFMD magnetic ordering. Carbon inter- 
stitial impurity in octahedral position is shown by red (dark) 
circle. 

Further, we have performed calculations for Fe32C su- 
percell with carbon in the octahedral void (Fig. [I} . This 
concentration is close to the eutectic point (3.6 at.%) at 
the phase diagram Fe-C 24j which is the most interesting 
from the point of view of metallurgy. Calculations of the 
total energy and magnetic moments, as well as optimiza- 
tion of positions of all atoms in the supercell have been 
carried out for FM, AFM, and AFMD structures (see 
Table Hi)) . Calculated exchange interactions for the first 
and second Fe-Fe bonds have shown in Fig. [2^,. These 
parameters agree well with previous calculations for fee 
iron [B| . Carbon in the octahedral void, even without rc- 



TABLE II: 

Lattice parameters, tetragonal deformations, magnetic mo- 
ments for nearest neighbors (nn) and next nearest neighbors 
(nnn), and total energy differences per iron atom for Fe32C 
unit cell; numbers in parenthesis calculated within LMTO 
values 





FM 


AFM 


AFMD 


a, A 


3.73 


3.56 


3.61 


c/a, nn 


0.94 


0.98 


0.93 


c/a, nnn 


0.99 


1.04 


1.00 


c/a, bulk 


0.97 


1.04 


1.00 


M, [ib nn 


2.3 (2.1) 


1.6 (1.8) 


1.9 (1.9) 


M, hb nnn 


2.8 (2.7) 


2.2 (2.4) 


2.6 (2.5) 


M, [i B bulk 


2.7 (2.5) 


1.9 (2.2) 


2.1 (2.3) 


E — Efm, meV 





47 (29) 


16 (13) 



laxation, changes the sign of nearest-neighbor exchange 
parameters from AFM (J x =-83 K) to FM (Ji=+96 K). 
Another effect is an essential increase of the next- nearest- 
neighbor exchange parameter (J 2 =48 K in Fig. [2k and 
J2=78 K in Fig. [2Jd). Interestingly, the relaxation makes 
the next-nearest-neighbor interactions even stronger than 
the nearest- neighbor ones (Fig. [2t). Probably, the gain 
of magnetic energy related to this effect is one of the driv- 
ing mechanism of the local tetragonal distortion. The 
main magnetic characteristics calculated in the SIESTA 
and in the LMTO are similar which confirms that our 
exchange parameters are reliable enough, at least, for 
qualitative discussions. 

It turned out that, in contrast with the case of pure 
7-iron, the FM ordering has the lowest energy in the pres- 
ence of carbon. The exchange parameters [18| calculated 
for the FM configuration presented in Fig. [5] also confirm 
that this magnetic configuration is stable. The accuracy 
of the Heisenberg model estimated from the difference of 
exchange parameters and values of magnetic moments in 
the FM and AFMD state is in the limit of 25% (see Table 
P. 

The mechanism of FM state stabilization by the car- 
bon impurity can be understood by investigation of Fe- 
C chemical bonding. We present the density of states 
(DOS) for Fe in Fe 32 C supercell, together with the local 
DOS for carbon impurity, in Fig. [3J One can see that 
for a broad energy interval (±2eV) near the Fermi level 
the hybridization of the sp-states of carbon with the d- 
states of iron is much more pronounced for the FM state 
(Fig. W) than for both AFM ones (Fig. [3p,c). This can 
lead to the energy stabilization of FM states in fee Fe- 
C alloys, which results in positive exchange interactions 
even without structure relaxation (Fig. [2p). The effect 
of anisotropic structural relaxation increases formation 
of strong FM bonds (Fig. [2p) and reduce the total iron 
DOS at the Fermi level (Fig. [3^) m comparison with the 
original AFM states. 
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FIG. 2: (color online) Exchange parameters (in K) for differ- 
ent Fe-Fe pairs in original fee lattice (a); in fee lattice with 
carbon interstitial impurity without (b) and with (c) relax- 
ation taken into account. Arrows indicate direction of atomic 
displacements during the relaxation. 



The solution enthalpy of carbon in 7-iron has been 
calculated from the total energies of FM Fe32C and of 
AFMD fee Fe (which have the lowest energies among trial 
magnetic configurations), together with the ground state 
energy of graphite. The result is 0.55 eV whereas ex- 
perimental value is about 0.4 eV @. Keeping in mind 
that ab initio calculations without taking into account 
local distortions and correct magnetic ground state give 
just a wrong sign for this quantity p} one can say that 
the agreement is rather good. Actually, this is even bet- 
ter since our calculations have been done for high enough 
carbon concentration and thus a fictitious carbon-carbon 
interaction presents. Estimations of this effect according 
to the standard elasticity theory (SBJ gives a value of or- 
der of 0.1 eV which should be subtracted from our result. 

We have done also calculations for the supercell FeiosC 
with the AFMD magnetic configuration. Starting from 
third coordination sphere of the carbon atom the tetrag- 
onal deformation c/a — 1 is approximately 7% which is 
close to the value for pure iron (see Table fl} . The cor- 
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FIG. 3: (color online) Total DOS of iron atoms (solid red 
line), DOS of iron atoms from the first shell (dotted blue 
line), and for carbon atom (dashed green line) for relaxed 
Fe32C supercell with carbon in octahedral void for different 
magnetic configurations. 



responding values of local tetragonal distortions for the 
first shell near the carbon impurity is 3% and for second 
shell is already 6%. This means that respective to the 
tetragonal deformation of the host carbon produces local 
distortion of its nearest surrounding at approximately - 
4% which is a bit smaller than for Fe32C supercell. The 
values of magnetic moments (in [is units) for the first 
three coordination spheres are 1.9, 2.5, and 2.6, respec- 
tively. Thus, local tetragonal distortion around carbon 
atom can be considered, indeed, as an effect of single 
impurity. 

Interactions between carbon atoms via lattice distor- 
tions (deformation interactions) determine decomposi- 
tion and carbon ordering processes in steel which are 
important for microstructure formation 241] . Octahedral 
voids in bec lattice are asymmetric themselves which re- 
sults in local tetragonal distortions around carbon in- 
terstitial impurity and rather strong deformation inter- 
actions. It is commonly accepted that the deformation 
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interactions in fee host are much weaker since the voids 
are symmetric and the interstitial carbons are consid- 
ered as purely dilatations centers (25|. We demonstrate 
that proper including the magnetic effects leads to the 
local tetragonal distortions around carbon in 7-Fe of the 
same order of magnitude as in a-Fe and, thus, traditional 
views on the importance of the deformation interactions 
in different phases should be reconsidered. 

These local deformations are intimately connected 
with the effect of carbon on local magnetic configura- 
tions of iron, ft turns out that carbon changes signs of 
some exchange integrals from AFM to FM. Similar effects 
of strong distortion dependence of exchange interactions 
has been discussed earlier for fee Fe-Ni alloys [5]. We 
show here that in addition to above mentioned distance 
dependence of effective exchange interaction the effect of 
Fe-C chemical bonding is also important. 

It should be noted, that our calculations are done for 
the ground-state case whereas an interesting temperature 
interval for Fe-C steels is above I0 3 K. Nevertheless, the 
investigated magnetic effects can be very important for 
understanding of structural distortions in 7-iron alloy. 
Indeed, local magnetic configurations and thus local dis- 
tortions can survive till relatively high temperatures. For 
the classical Heisenberg model on fee lattice, the mean- 
field estimation for the Curie temperature is equal to 
the energy difference between FM and AFM configura- 
tions which gives a value of order of 500 K (see Table 
HI)) . Quantum effects for magnetic moments of order of 
2 fiB increases this estimation by a factor of 2 18, [26| 
which allows to assume that, at least, up to 10 3 K local 
magnetic correlations will survive. This temperature can 
be higher for higher concentration of carbon. There are 
some direct experimental evidences that magnetic effects 
are important for the austenite-to-ferrite (a — 7) trans- 
formation in steel [27j |. 

In conclusion, the complex magnetic state with strong 
tetragonal distortions is predicted for 7-iron near carbon 
impurities. The calculated exchange interactions show 
the strong tendency to formation of local FM clusters. 
This effect changes drastically carbon-carbon deforma- 
tion interactions in 7-phase and thus should be relevant 
for the martensitic transformations in steel. 
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